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In this paper the structure of the next-to-leading (NLO) static gluon self energy for an anisotropic 
plasma is investigated in the limit of a small momentum space anisotropy. Using the Ward identities 
for the static hard-loop (HL) gluon polarization tensor and the (nontrivial) static HL vertices, we 
derive a comparatively compact form for the complete NLO correction to the structure function 
containing the space-like pole associated with magnetic instabilities. On the basis of a calculation 
without HL vertices, it has been conjectured that the imaginary part of this structure function 
is nonzero, rendering the space-like poles integrable. We show that there are both positive and 
negative contributions when HL vertices are included, highlighting the necessity of a complete 
numerical evaluation, for which the present work provides the basis. 
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I. INTRODUCTION 



The difficulties in explaining the strong collective phenomena observed at the Relativistic Heavy Ion collider (RHIC) 
[H by perturbative QCD at finite temperature have led to extensive studies of the consequences of the inevitable 
presence of non-Abelian plasma instabilities 0, 0, [3| in a plasma with momentum-space anisotropy. To leading order 
in the coupling and for small gauge field amplitudes, the dynamics of plasma instabilities are determined by the 
generaHzation of the hard-thermal- loop (HTL) [1,0,01 gauge boson self-energy to anisotropic situations [1, H, [13, [lH . 

In equilibrium, it is well known that non-bilinear terms in the HTL effective action [H, [3| are important at 
next-to- leading order (NLO). In order to obtain complete NLO corrections to HTL dispersion laws such as plasmon 
damping constants, one must resum both propagators and vertices fisljl^. [itI. fT§. [l9| . In anisotropic systems, the 
non-bilinear terms in the corresponding hard-loop (HL) effective action [20| are important for the dynamics of non- 
AbeHan plasma instabilities at large gauge field amplitudes, which has been studied using real-time lattice simulations 
[2ll . m, m, 0; m, HEI- up to now, analytic calculations at NLO in anisotropic systems have been considered 
prohibitively difficult. In addition to the obvious technical difficulties associated in dealing with expressions that 
contain huge numbers of terms, these calculations contain new conceptual problems arising from the fact that the 
anisotropic gluon propagator contains non-integrable space-Hke poles. 

In perturbative estimates of jet quenchin g an d momentum broadening in the anisotropic quark-gluon plasma, such 
non-integrable singularities also appear (27l. l28| . It has been suggested in [131 that the static gluon self energy might 
develop a non-zero imaginary part at NLO that regulates these singularities. In thermal equihbrium, the imaginary 
par t of the static gluon self energy, which is an odd function of the frequency, vanishes due to the KMS condition 
|29| . In the anisotropic case, it seems possible that there is a finite, discontinuous contribution. Ref. [23I provided 
support for this conjecture in the form of a partial calculation of the anisotropic NLO static gluon self-energy. This 
calculation included however only a tadpole diagram with a bare four-vertex. 

In this paper, we provide the basis for a complete analytic calculation of the anisotropic NLO static gluon self- 
energy. In the HTL case, only static loop momenta have to be considered, no fermion loops contribute, and all HTL 
vertex corrections vanish js^l- However, the anisotropic HL vertex corrections do not vanish in the static limit [2o| . 
The relevant diagrams are shown in Fig. [T] (since the ghost self energy vanishes at leading order, the ghosts do not 
need to be resummed). The solid dots indicate leading order propagators and vertices which are obtained from the 
hard loop (HL) effective action [23|. The third diagram is the HL counterterm which must be subtracted to avoid 
double counting. 

In this paper, we give an analytic result for the integrand corresponding to the diagrams in Fig. [TJ After extensive 
algebraic manipulations, the final expression has a relatively simple form, because of cancellations that occur between 
different contributions from the tadpole and bubble contributions. We divide the result into 'tadpole-like' contributions 
(which contain only one HL propagator) and 'bubble-Hke' contributions (which contain two HL propagators). We 
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calculate the contribution from 'tadpole-like' contributions with a bare vertex. We compare our result with the 
(corrected) result of [IS] for the tadpole diagram (without taking into account cancellations with the bubble diagram) 
with a bare vertex. 




FIG. 1: The diagrams that contribute to static gluon self energy. All lines correspond to gluon propagators. The dots indicate 
hard loop propagators and vertices. The cross denotes the HL counterterm. 



II. NOTATION AND HARD ANISOTROPIC LOOPS 



In this section we define our notation for the anisotropic HL quantities and give the explicit form of the static HL 
propagator in Feynman gauge. 

At zero temperature, field theory can be formulated covariantly. At finite temperature, covariance is broken by the 
vector — {1, 0, 0, 0) which specifies the rest frame of the thermal system. For anisotropic systems we need (in the 
simplest case) one additional vector to specify the direction of the anisotropy. We consider the case in which there is 
one preferred spatial direction along which the system is anisotropic (in planes transverse to this vector the system is 
isotropic). In the context of heavy ion collisions, we can take this direction to be the beam axis (z) along which the 
initial expansion occurs. 



A. Distribution Functions 



We define the isotropic particle distribution: 

fisoip) = 2Nf[nip) + n{p)] + AN^n^ip) (1) 

In equilibrium, we have 

We define the Debye mass from the equilibrium distribution^ 



2 2 

rriD = g 



Following [1, we can construct an anisotropic distribution from any isotropic distribution of the form /iso(p^) 
by writing 

f{p) = fisoip^ +ap-sf) (4) 

where ^ > —1 is the anisotropy parameter. A value ^ > corresponds to a contraction of the distribution and 
> ^ > — 1 corresponds to a stretching of the distribution. In this paper we restrict ourselves to weakly anisotropic 
systems for which |^| <C 1. 



[1] In Ref. @, HOl Eq. GJ differs by a factor of 2 and Eq. |[TJ differs by a factor 1/2. The definition of the Debye mass is the same. 
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B. Projection Operators 

Following [3l| we construct the vector 

n\k) = {6'^ - k'y/k^)5^^ = S'^ ~ Ffc(3) ; fc(3) fc . z = k^^^ (5) 

which satisfies n(fc)* /c* = 0. Throughout this paper we will frequently use the indices {fc, q, r} to denote momentum 
arguments. For example: we will write n\ :— n'^{k). We also use Latin letters l,t,- ■ ■} to denote spatial indices, 
with the exception that the indices fc, q, r are reserved and used exclusively to denote momenta. Using the vector in 
([5]), we construct the projection operators: 

pl_c__hkj_. 2_M2. p3 _ _Ml _ 4_ A. . P4 - + ■ p5 - pi ^ ((\\ 

which satisfy: 

PiiPij — Pij ; PiiPij — ; PaPfj — Pfj ; PaP^j + PaPij = Pij (7) 

p2 p2 p2 . p2 p3 fi . p2 p4 , p4 p2 p4 

p3 p3 p3 . p3 p4 ft 

p2 p5 n . p3 p5 n . p4 p5 p4 . p5 p5 p5 

^iPij ~ ^iPij ~ kikjP^j = kiP^j = 

n^i^i^n^; - n^/>^ = ; n^n^/^=0; = 

Tr = 2 ; Tr = Tr = Tr = 1 ; Tr P^ = 

C. Lowest Order Self Energy 

We define the polarization tensor using the relation: 

z?;J(i^) = iD%y' - n^. = -{g,..K' + n^.) (8) 

The HL gluon self energy is gauge invariant and satisfies the usual Ward identity: K^^W^u — 0. Consequently, we 
only need to calculate the spatial components. At finite temperature, there are two independent components which 
are called the transverse and longitudinal parts. For anisotropic systems the self energy can be decomposed into four 
independent structure functions: 

n^. = 4 au + Pi h + P^ Ik + PtjSk (9) 
The structure functions are calculated from: 

X = P(2:)y ; xe{a, ^, 7, 6} (10) 

P(a),, =Pg; rm.^P^; Vi-f)., ^ ^ 2P^ ; V{5) ^ P^ / {2enl) 
In the limit of small anisotropy parameter ^, the lowest order results are (9I. [lo|: 

^Wd , c ( '^''k) kWo 1.-, 2n 2 V . ( {2 ~ 3nl) TT^koml Trkoml\ 



fc2 ^ I 3fc2 3' I \ 8fc 4fc 

{3nl - 1) Cfcg^D H^l 



3fc2 fc2 



7fc = M + o^D^k 



3fc2 ' 3 ^ V 4fc 



3/ 

For convenience we define 



^'^-^ 3F + 4P (^^^ 



fco4 := 4 ; k^fik ■■= h (12) 
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D. Propagator 

We calculate the static propagator in the covariant gauge (with Feynman gauge parameter) by inverting ([8]). This 
inversion is complicated by the fact that the anisotropic propagator depends on the two fixed vectors (1,0,0,0) and 
(0,0,0,1), as well as the 4-momentum K^. The method is described in [3l|. The result is: 

A, (fc) - ~P?, Y2 ~ ^■4(fc)Pf, + nlFu {k)P^, (13) 

Doi{k) = Fsinkji 
Doo = Fcik) 

FB{k)=6kGB{k) 

k'Fcik) = +ak+ik)GB{k) 

nlFoik) = PkCeik) 

G-j^\k) = e + ak 

Gs\k) = (1 - Pk){k^ + ak + 7fc) + k^nlSl 

In the small-^ limit, the result (11) implies that Gb has a pole at negative k^ — —rnjj + 0{^), which corresponds 
to the usual electric (Debye) screening. For non-vanishing ^ space-like poles also appear. Writing k ■ z — cosO, and 
n\ — l~(k-zY = sin^ 9, we find that these space-like poles appear at: 

GB(fc) : fcV™l5 = ^ai-2n2) = i^cos20 (14) 

GA{k) : e/ml = 1 ^ (1 - n^.) = i ^ cos^ 9 

When ^ is positive, corresponding to an oblate particle distribution, GB{k) has a space-like pole for tt > 6 > 3tt/4 
and 7r/4 > 9 > 0. For ^ negative, the space-Hke pole occurs for 37r/4 > 6 > 7r/4. GA{k) has a space-like pole for any 
positive value of ^ unless fc • z = 0. These poles of Ga correspond to the magnetic Weibel instability j32] . 



E. Vertices 



In this section we give our notation for the HL vertices [20|- We define: 



9^ f df 



2 Jp dPl^ 



(15) 



The 2-point function is: 



= .9 



up 



PQ 



(16) 



The 3-point function is 



;„2 



:= 



P- K PQ PQ PR 



(17) 



[2] In Ref. Q there is a missing factor (-1) in Eq. (35) 



5 



The 4-point function is 



Kbcf - -2g^{Xabcd M^"^" + 2 cyclic permutations) ; Xatcd = tr (T4T„ [n, TJ]) (18) 



P ■ K P -Q P ■ {Q + L) P -Q P ■ L P ■ {L + S) P ■ L P ■ S P ■ [K + S) 

Momenta are taken to be incoming. When the momentum arguments are not written exphcitly, they are taken to be 
in the order {K, Q,R = —K — Q) for the 3-point function and {K, Q,L,S = —K — Q — L) for the 4-point function. 
The HL vertices satisfy the Ward identities: 

if^r^'^^(if, Q, R) = n'-^Q) - n"^(i?) (19) 

K^M^'^^'^iK, Q, L, S) = V^^iQ, L, -L-Q)- ^"^"{-1 - S, L, S) 
The tadpole form of the 4-point vertex has a simpler form: 

KbcciQ^ -Q^ K, -K) :- 2.92 CASab M'^-'^^Q, -Q, K, -K) (20) 
M^^'^'^iQ, -Q, K, ~K) = -iipP^'P^P^P" 



P-K P-Q{{P-Kf-{P-Qf) j 
This vertex satisfies the Ward identities: 

KxM>"'^''{Q., -Q, K, -~K) ^ ^2T'""'{K, Q, ~K - Q) (21) 
RxK^M^-'^'^iQ, -Q, K, -K) = 2W(-K - Q) - 2W {Q) 

We also define the bare vertices: 

(ro)r' = *5./afccrr (22) 

{MoTabcciQ^ -Q^ K, -K) = 252 CASab M'^'^^^Q, -Q, K, -K) 
M^"^" = -g^^g"" - g^^g""" + 2g^'".g^'' 

Using these definitions we write the integrand corresponding to the first two diagrams in Fig. [T] as: 

1 f d^k II 

n^"(Q)= 2 5'CA<5,bT j j^[{T'' + Tf^^-Dxy{K){T'' +Tf''^ Drr'{K) + {M^ + MY-^"Dxa{K)\. (23) 

The order of the momentum variables is {K,Q,R) for the 3-point functions and {Q,~Q,K,~K) for the 4-point 
functions. Eq. l(23|) is to be understood as restricted to the static case and small momentum anisotropy such that the 
prefactor of the integral is simply the temperature T of the isotropic distribution function ([2]) . 



III. INTEGRAND FOR ^nlo 

We calculate the NLO contribution to a which is obtained from lfTO|) as: 

aNLO =i'fj(<Z)%(q), (24) 

with Ily (g) as given in l(23|) . The resulting integrand can be rewritten in a comparatively compact form. The method 
is similar to the calculations of the integrands that give the NLO equilibrium gauge boson polarization [l3|, and 
the NLO equilibrium fermion self energy [15, 16, 18]. The basic procedure is as follows. We substitute the static 
propagator as given in Eq. lfT3|) . Then we use the Ward identities in Eqs. lfT9|) and (pTI) . This produces factors 



[3] In Ref. |20l | there is a missing factor (-1/2) in Eq. (39) and a missing i in Eq. (38) 
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proportional to components of the HL polarization tensor in the numerator. These factors can be expressed as pieces 
of inverse propagators, which cancel with the corresponding components of the original propagators. This procedure 
allows us to identify cancellations between various pieces of the first two diagrams shown in Fig. [TJ We divide the 
result into several types of contributions: (i) terms in which all factors of the HL propagator have cancelled; (ii) 
tadpole-like terms which contain only one propagator; (iii) bubble-like terms that contain two propagators and no 
HL vertices; (iv) bubble-like terms that contain two propagators and one HL vertex; and (v) bubble-Hke terms that 
contain two propagators and two HL vertices. We thus write 

aNLO = ^ 5' <5., T I ^ [ g ] (25) 

where the result for the various contributions to the integrand /^*^ are given below. To compactify the expression we 
use the notation: 

VkJ{k,q):=f{k,q) + f{q,k) (26) 

Hkqrfik, q, r) := /(fc, q, r) + /(g, r, k) + /(r, k, q) 

F{k,nk,fJ.,iy) = 9^39.3FD{k) - {9^39.0 + g^og.s) Fsik) +9^o9^oFc{k) 

The separate integrands then read: 

(i) terms in which all factors of the HL propagator have cancelled: 

,(1) _ _ 2(g^+^-fc) ^ 2(2r ■ fc ^ Sfcara) 

^ ~ ;„2„2 + ;.2„2„2 y^^') 



(ii) tadpole-like terms: 



An (h\'u ( g-fc , jq-k-k^qs) q^r^ 

\ q'^r'^ q^r'^n^ q'^r'^nj,/ 

1„ /M3333 . 2Mii33 



GA{k)Vkq [^^ + M, 



2--v^.- ^ ^2^2 "JJ „2 

^ F (fc, nfc, v) ((M,„33 - M,„iin2) q^) 



9 9 
q'^n'^ 



-q 

(iii) bubble-like terms with no HL vertices: 

-GA{k)Vkq^k3q3&k{f^qFc(r)~2^^FB{r)) (29) 



2 

-G A{kyPkq-r^k3q3^klqFD(T) 



+GA{k)Vkq\ (nl {nl + 1) - 2nl) (q^S^Fsir) - iqFu(r)) 

"^k 

+^q^nl5qFB{r) {Fc(k) - Foik)) 
+V/3^ (FB{k)FB{r) ~ Fc{k)Fc{r)) 
+Anl^q {FD{k)FD{r) - FB{k)FB{r)) 

+GA{k)GAir)n, r + ^^'^^ + 2 jr' {nj - l) - g^n^) _ 2k^\ ^^^^ 

I /j2^2 J^2j.2j^2^2 J^2j.2^2 ^2^,2^2 I 
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+GAik)GA{r)'Jlkqr^ {{nl - l) - nl {nl - l)) 7, 



(g2 + 8r • fc n2 - 4A:3r3) Fc(A:)Fc.(r) 



GA(fc)FB(r)7'fc, 2 1 - ^ + 10 1 - — {k'ni - r'n^ 



+GA{k)GA{r)'Rkqr 5k'ni 



(iv) bubble-like terms with one HL vertex: 



1 1 



2q^ 



= GA{k)GAir)nkgr^^q3 (r333 " n^TrFg) 



(fc, nk,fi, A) F (r, n^, i^, t) {g^og^o - g^zg^zn^) (fcs - fs) r;,,3 

+GA{k)Vkq\'F (r, rir, 3, A) ((93 - ^-a) Fass - n^^qzivTx) 



(v) bubble-like terms with two HL vertices: 



(30) 



= GA{k)GA{r)nkgr [- 



^ 333 



1 F3 • F3 



2r„3 • F 



3 ■ i ,^3 



-GA(fc)F (r, n^, v) Vkq 



+¥ (fc, rifc, ^, v) F (r, n^. A, r) F^a ■ F 



F • F 



F„.3F 



F33 • F33 

F,i33F„33 

nlnl 



(31) 



IV. IMAGINARY PART OF TADPOLE-LIKE TERMS 



We have evaluated numerically the imaginary part of the tadpole-Hke terms given in l(28|) with the HL vertices set to 
zero. The result is given in Fig. [2l We compare with the result of [131 for the tadpole diagram with a HL propagator 
and a bare vertex (see Appendix A): 



V3 ■ 

Im aNT.o = 

" 10247r-^ 



■Ty^niDNc (56 + 5\/21og (3 - 2V2^ - 5\/21og (3 + 2Vtj^ 



(32) 




FIG. 2: (Color online) The result for onlo as a function oi qs/q for the simple tadpole (the (red) straight line) as given in Eq. 
I|32p . and for the 'fuU'-tadpole (the (blue) curve) as given in Eq. dSHl). The factor N^^mDg'^T/{16Tv'^) has been extracted. 
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By considering only the tadpole diagram with a HL propagator but a bare vertex, Ref. [27| has obtained an 
imaginary part Im a which is of the order gm|) and strictly negative. Our result, which includes also contributions 
from the bubble diagram with HL vertices where the latter have cancelled one of the two HL propagators, shows that 
this finding is not generic as far as the sign of the result is concerned. Evidently, the complete static NLO gauge 
boson self energy contains positive as well as negative contributions to the imaginary part, and it remains to be seen 
whether a complete numerical evaluation of the expressions that we have presented here will lead to a nonzero result. 

V. CONCLUSIONS 

The main result of this paper is a relatively simple result for the integrand that gives the complete NLO contribution 
to the a component of the gluon self energy, in the static limit. We have calculated analytically the contribution to the 
imaginary part from 'tadpole-Hke' terms, by which we mean contributions which involve only a single HL propagator 
after Ward identities have been used to reduce the diagram involving two HL propagators and two HL 3- vertices. The 
result is of the order predicted in [23|, but the tadpole-like contribution is not strictly negative. 

Indeed, from previous experience with equilibrium hard thermal loop calculations, one would expect cancellations 
between the various contributions to the full integral. For the equihbrium case, the KMS conditions actually guarantee 
that the imaginary part of the self energy is zero in the static limit, at all orders. Out of equihbrium, the result must 
be odd in the frequency, which means that any nonzero result must contain a discontinuity at vanishing frequency. 
The existence of this discontinuity at the order grr?j~, has been conjectured in |23| and used for a perturbative estimate 
of the order of magnitude of jet quenching and momentum broadening parameters [13, US]- In order to determine 
whether there is in fact a non-zero contribution to the imaginary part of the complete static gluon self-energy at 
next-to-leading, the full integral, as given in Eqs. l(28|) - (|3T|) . must be evaluated numerically. This calculation, which 
involves some nested integrations with highly nontrivial integrands, is in progress. 

APPENDIX A: COMPARISON WITH REF. [131 

In this Appendix we point out a few errors that occurred in the calculation of the tadpole diagram with bare vertex 
which was presented in the Appendix of Ref. [23|. The main mistake is the incorrect assumption of a relation of 
the form J'^H*-' c>c n^n^n*-' which leads to the incorrect simplification: 

/? = (Al) 
2a + 7 = %ny. (A2) 

However, the results given in Ref. [131 do not satisfy (|A2p because of typographical errors. We write below the results 
from [i^l for onlo and 7nlo from Eq. (A3) of that paper, with an extra factor -\/2 in front of the logs, and an extra 
factor (-1) in front of the result for 7nlo- These extra factors are written in square brackets. 

«NLo[2l = "I^^'^^"'^^^ " (3 + 2V2)) {2 I) (A3) 

I^7^Lo[2l = ["^^WIV^'^^"''^^'' (2- [V2]ilog(3 + 2^/2)) 
These 'adjusted' expressions satisfy Eq. HA2p . However Eq. (|Aip is not valid (not even at leading order). 
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